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Quantum gases of atoms and exciton-polaritons are nowadays a well established theoretical and 
experimental tool for fundamental studies of quantum many-body physics and suggest promising 
applications to quantum computing. Given their technological complexity, it is of paramount interest 
to devise other systems where such quantum many-body physics can be investigated at a lesser 
technological expense. Here we examine a relatively well-known system of laser light propagating 
through thermo-optical defocusing media: based on a hydrodynamical description of light as a 
quantum fluid of interacting photons, we investigate such systems as a valid, room temperature 
alternative to atomic or exciton-polariton condensates for studies of many-body physics. First, we 
show that by using a technique traditionally used in oceanography, it is possible to perform a direct 
measurement of the single-particle part of the dispersion relation of the elementary excitations on 
top of the photon fluid and to detect its global flow. Then, using a pump-and-probe set-up, we 
investigate the collective nature of low-wavevector sound modes of the fluid and observe signatures 
of superfluid behaviour. 


INTRODUCTION 

Quantum gases, i.e. systems in which the thermal 
de Broglie wavelength is larger than the average parti¬ 
cle distance are an ever increasingly important area of 
theoretical and experimental study with applications as 
diverse as quantum computing and quantum simulation 
of general relativity models. The most important ex¬ 
ample of such physics is the Bose-Einstein Condensate 
where a macroscopic number of bosonic particles at low 
temperature share the same wavefunction. After early 
work on liquid Helium, the physics of condensates has 
been experimentally studied in ultracold atomic gases 
and, more recently, in exciton-polariton fluids in semi¬ 
conductor microcavities. Building on these latter studies, 
an ever growing community is now active on the study 
of the so-called quantum fluids of light, where the many- 
photons forming the beam are seen as a gas of interacting 
particles via the optical nonlinearity of the medium [T]. 
Among the many hydrodynamic effects that are presently 
being studied in such fluids of light, we can mention tur¬ 
bulence [ 2 ] , where the general physical processes found in 
a fluid or superfluid can be found in many other systems 
such as plasmas or in astrophysical systems. Another in¬ 
triguing direction is so called ’’analogue gravity” where 
flowing fluids are used to model gravitational spacetime 
geometries nn. Waves propagating on top of a flowing 
medium behave the same way as light waves in a grav¬ 
itational field. Hence, by tailoring the properties and 
flow geometry of such a medium, it is possible to mimic 
gravitational spacetimes in a laboratory experiment. As 
an alternative to actual flowing fluids, optical analogues 
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obtained using laser beams propagating through nonlin¬ 
ear self-defocusing media or confined in nonlinear opti¬ 
cal cavities have been proposed to create a photon fluid 
with the desired flow pattern. In these systems the fluid 
properties can be controlled through the phase and in¬ 
tensity of the incident optical field as well as by the re¬ 
fractive index profile and the structural boundaries of 
the medium [zHg. For example, a theoretical framework 
for photon fluids [n H [g cni cn has shown that sound¬ 
like excitations experience a curved spacetime and hence 
are promising models for analogue gravity experiments. 
A full quantum simulation showing evidence for analog 
Hawking radiation from a black hole horizon in a fluid of 
light was performed in m and the first experimental re¬ 
ports of black hole horizons in fluids of light were reported 
in Oils]. The effective photon-photon interaction in the 
nonlinear medium arises from the third-order Kerr non¬ 
linearity [IS],i.e. the local change in refractive index pro¬ 
portional to the light intensity An = n2\E\‘^: a defocus¬ 
ing nonlinearity with 712 < 0 giving repulsive interactions 
is required to observe a dynamically stable photon fluid. 
Such a nonlinearity may be obtained through the thermo¬ 
optic effect, i.e. an intense laser beam locally heats the 
medium, which reacts by decreasing the refractive index 
proportionally to the deposited laser energy thus giving 
rise to an effective repulsive photon-photon interaction 
m- Here, we study and propose a relatively well known 
nonlinear medium, i.e. the thermal defocusing medium 
described above however, revisited within the context of 
a fluid of light. The aim is to study the extent to which 
such systems can indeed represent a novel platform for 
studying, at room temperature, the properties of quan¬ 
tum gases with an emphasis on (but not limited to) the 
exciting perspective of reproducing superfluid behaviour, 
phonon propagation according to the Bogoliubov disper¬ 
sion, and particle generation from a black hole horizon 
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and in an analogue gravity set-up. In the propagating ge¬ 
ometry considered here, the dynamics of the fluid takes 
place in the transverse plane (x,y) of a laser beam prop¬ 
agating through a bulk nonlinear medium, so that the 
propagation coordinate z plays the role of the time co¬ 
ordinate t [TOl [TT] . Hence, by imaging the evolution of 
the beam profile along the z-coordinate, it is possible to 
detect the spatio-temporal evolution of effective sound 
modes in the photon fluid (i.e. small oscillations on the 
transverse beam profile) and extract information on the 
dispersion relation of its elementary excitations. Using 
a technique inspired from oceanography, we obtain in¬ 
formation on the high wavevector part of the dispersion 
where excitations have a single-particle nature and we de¬ 
tect the effect of global flow of the fluid. Then, a pump- 
and-probe technique is presented, where the probe beam 
is used to create elementary excitations in the fluid gener¬ 
ated by the pump. Information on the elementary excita¬ 
tion dispersion is obtained from the lateral shift of the in¬ 
terference pattern of pump and probe beams: signatures 
of the collective nature of the excitations and of the non¬ 
locality of the photon-photon interactions are identihed 
and discussed. An important result of this study is the 
experimental verification of the existence of a (phonon) 
wavelength regime in which, using the terminology in¬ 
troduced by Chiao et ah, photon superfluidity can be 
observed [SI- 


THEORETICAL MODEL 


The propagation of a monochromatic laser beam with a 
vacuum optical wavelength A in a self-defocusing medium 
can be described at steady state by the non-linear 
Schrodinger equation (NLSE): 


Zk no 


(1) 


Where z is the propagation direction, k = 27rno/A the 
optical wave vector and tiq is the linear refractive in¬ 
dex of the medium. The Laplacian term with re¬ 

spect to the transverse coordinates Y{x,y) describes the 
linear diffraction of the laser beam, where the second, 
nonlinear, term proportional to 712 refers to the optical 
Kerr nonlinearity, i.e. an intensity-dependent refractive 
index. In this work, we will always take 712 < 0, that is a 
self-defocusing nonlinearity, which guarantees transverse 
stability of the beams. By writing the held E in terms of 
its amplitude and phase E = Ebgexp{i4>), equation Q 
can be re-written as two differential equations |12j 


drP + '^{pv) = 0 
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= 0 
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formally identical to the density and phase equations of a 
two-dimensional BEC, the temporal variable r = zn^jc 


being proportional to the propagation distance. The op¬ 
tical background intensity is identihed as the pho¬ 

ton huid density p and the phase </> dehnes a huid velocity 

V = {c/kno)^(j) = Vip. 

The small amplitude perturbations E = E^g + e (E^g 
is the background held) can be described within the Bo- 
goliubov theory in terms of sound waves on top of the 
photon huid [101 [IH]. In the spatially homogeneous case 
where E^g does not depend on the transverse coordinate 
r, a sound mode of wavevector K and amplitude aK has 
a plane-wave form 


e = aKUK exp {—iflr + iKr) -\- vk exp (iflr — iKr) 

(4) 

and its angular frequency H (in the temporal r variable) 
satishes the dispersion relation 






(5) 


Here, c is the speed of light and v the background 
how velocity. Using the terminology of hydrodynamics, 
we note that for low frequencies the dispersion of the 
sound modes has a sonic character H oc K, whereas for 
higher frequencies the second term dominates and the 
sound modes follow a single-particle quadratic dispersion 
U cx K^. The separation between these two regimes 
dehnes a characteristic length ^ = A/2-^7io|772| 
usually called the healing length in the BEC litera¬ 
ture m- As a result, only long-wavelength sound 
modes with A ^ ^ follow a sonic dispersio n with a 
constant sound speed Cg = ^/c^\n2\JEbg\^7^■ This is 
of particular relevance, since only a linear dispersion 
guarantees superhuid behaviour |10j and, in the context 
of analogue gravity, Lorentz invariance in the acoustic 
metric 13 uni- 


Translating these concepts in to the language of op¬ 
tics [To] , K = is the magnitude of the trans¬ 

verse wave vector of the light held and H = {c/no)AKz 
is related to the change Ak^ in the wave vector along 
the propagation axis. To experimentally access the low 
frequency modes i.e. the A ^ regime, the light has to 
propagate in the nonlinear medium for at least one oscil¬ 
lation period T, so with A = CgT ^ one arrives at a 
minimum propagation length given by Zmin ~ cT/uq ^ 
X/An. We therefore see that for a reasonable nonlinear¬ 
ity of An = 10“'*, long samples with at least 10 — 20 cm 
length are required to observe the propagation of modes 
that lie in the linear regime of the dispersion. 


NON-LOCAL MEDIA 

A natural choice for the nonlinear material is therefore 
a liquid that exhibits a thermo-optic nonlinearity. The 
main advantage is that samples of the required dimen¬ 
sions can easily be built and that they are isotropic. In 
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thermo-optical media, self-defocusing of the laser beam 
arises due to partial absorption of the beam, which re¬ 
sults in heating effects and thus to a decrease in the re¬ 
fractive index nnun]. Thermal nonlinearity is usually 
highly nonlocal, i.e. the change in refractive index at any 
position depends not only on the local intensity, but also 
on surrounding field intensities [HI [H] ■ This is a result 
of heat conduction inside the nonlinear medium and has 
a striking influence on the wave dynamics in nonlinear 
optics [24H26] . The medium response may be described 
by a response function R{r, z) so that An can be written 
as 

An{r,z)=j j R{r — y')I{v') dr'dz , (6) 

with 7 being a normalisation factor and / the field 
intensity. In general, R depends on the nonlocal process 
in the material and can be of different shapes. It is also 
important to note that due to the long-range nature of 
thermal diffusion, R will also depend on the boundary 
conditions, e.g. on the shape and materials used for 
the cell containing the liquid [23] ■ For thermo-optical 
media an isotropic exponential decay was proposed 
for R{r,z) = 1/(20-^) exp (— IcrL)^ where the 
width of the response function defines the degree of 
nonlocality [23IE2I. 

Straightforward extension of the Bogoliubov theory to 
non-local media leads to a modified Bogoliubov disper¬ 
sion of the form [niET]: 

R[K,n,n/c)K^+^K^ (7) 

where R is the three-dimensional Fourier transform of 
the response function R{r,z). It is easy to see that 
well within the paraxial approximation Kjk,vjc <SC 1, 
the rescaled Bogoliubov frequency n^Vl/c ^ K, so 
the main contribution to the non-locality comes 
from the K dependence of i?, and we can simplify 
R{K, n^n/c) ~ R{K, 0) = 1/(1 -k K'^al) in Eq. Q. 

Examples of the Bogoliubov dispersion in non-local 
media with different parameters are shown in Fig. [^ For 
very low wavevectors K 1/^, l/a^ the nonlinearity is 
responsible for the linear, sonic shape of the dispersion 
n ~ CsK. At very large wavevectors K 3> 1/^, I/ctl, 
the dispersion recovers the quadratic single-particle 
dispersion fl oc . The effect of the non-locality is most 
interesting for ctl 3> in this case, there is a wide range 
of I/ctl < K < 1/^ where the nonlocality suppresses 
the effectiveness of the nonlinearity and the dispersion 
is pushed towards the single-particle dispersion already 
at iF ~ I/ctl ^ l/'f. For a sufficiently small aL and a 
sufficiently fast decrease of R (faster than the one we 
are considering here, e.g. Gaussian), a ring of roton-like 
minima appears in the dispersion for finite K HZ). As a 
sonic shape of the dispersion of excitations is essential 



K (10® m'^) 

FIG. 1. Bogoliubov dispersion relation according to Eq.([^ 
and[^. Local nonlinear media (An = 10“"^) show a linear 
part for wavelengths A ^ and |An| > 0 (black solid line, 
^ ~ 350 /rm, equivalent to K= 2.7 x 10® 1/m). In nonlo¬ 
cal media the effectiveness of the nonlinearity is suppressed 
for wavelengths A < (t_l that results in a parabolic disper¬ 
sion relation (the dotted blue lines correspond to a nonlocal 
nonlinear dispersion with cri, = 1, 5 and 10 /rm). 


for gravitational or superfluid analogies to be rigorous, 
experimental studies in this direction require having 
control over the nonlinear and nonlocal properties of the 
material. 

In the following, we report experimental evidence of 
a photon fluid in a nonlocal, self-defocusing medium. 
We adapted a method from oceanography to measure 
the spatiotemporal dispersion relation and show that by 
controlling the phase and intensity of the laser beam, we 
can imprint a flow as well as tune the velocity of acoustic 
waves. We then demonstrate a second method, based on 
tracking the sound speed of waves in the photon fluid 
as a function of input power, that allows us to measure 
both the local nonlinear refractive index change. An and 
the nonlocal length of the medium, 


EXPERIMENTAL RESULTS I: 

OCEANOGRAPHIC TECHNIQUE 

The experiments were carried out by launching a 
collimated CW-laser beam with 532 nm wavelength 
through a cylindrical tube with length L = 13 cm (2 cm 
diameter) filled with a methanol/graphene solution. The 
beam is magnified to a waist diameter of 1.6 cm by a 4-f 
telescope and the central portion with a radius of 1 cm 
is selected with an aperture to ensure a relatively flat 
intensity profile. Methanol is known to have a negative 
thermo-optic coefficient of dn/dT « —4 x 10“^ 1/K 
but has a very low absorption coefficient of a = 0.0006 
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cm~^ [2S]. Nanometric graphene flakes (7 nm average 
size) are therefore dissolved in the medium in order to 
increase absorption. Differently from widely used dyes, 
graphene does not exhibit any efficient fluorescence 
mechanism so that most of the absorbed laser energy 
is converted directly into heat. By using a very dilute 
solution of graphene, it was possible to prepare samples 
with a total absorption that could be controllably 
chosen in the full 0 to ~ 100% region over the sample 
length and adjusted in our experiments so as to have 
20% absorption. A camera/lens system is mounted on 
a computer-controlled, longitudinal translation stage, 
so that the transverse beam profile can be imaged at 
different positions along the propagation axis of the 
beam. The imaging was performed at a constant image 
distance to ensure uniform scaling in transverse and 
z-directions during scanning. Since the z-axis maps 
into a time coordinate, it is possible to measure the 
spatiotemporal evolution of small intensity fluctuations 
on top of the beam by scanning along z and measuring 
the 2D beam profile at defined increments of Az. The 
input beam diameter was large enough so that linear 
diffraction is negligible along the sample. To directly 
measure the dispersion relation, a technique adapted 
from oceanographic studies was used where surface 
waves that naturally occur were recorded in both space 
and time and by calculating the Fourier transform one 
obtains the dispersion of the medium . 

In our experiment, the set of images are stored 
in a 3D (x,y,z) dataset: the origin of the observed 
fluctuations is due to (static) classical amplitude noise 
that is always present on the transverse profile of 
any laser beam (unless one invests significant effort 
to completely spatially filter it out). The dispersion 
AKz) can be calculated via the Fourier transform 
of a lineout of the measured intensity signal /(x, 0, z). 
The resulting spectrum can be converted to frequency 
using D = (c/no)AKz. To enhance the signal-to-noise 
ratio, the data was averaged over 200 line-outs along 
the y direction. Examples of the measured dispersion 
of the photon fluid are shown in Figs. |^a), (b) and (c). 
The spectra show some “folding” determined by the 
lattice of points at which the signal is sampled and then 
Fourier-transformed. The large signal at (K = 0, D = 0) 
corresponds to the background pump beam. Some of 
the spectra also show vertical stripes or other features 
that do not lie on the dispersion curves and which we 
were able to attribute to small defects present e.g. on 
the input or output window facets. The main finding of 
these measurements is the observation of an apparently 
purely parabolic dispersion for all spatial K-vectors and 
for all sets of parameters. Three different measurements 
are shown in Fig. for three different angles of the input 
beam with respect to the sample (and imaging) axis in 
the X — z plane. This tilt angle creates a linear phase 
gradient along the tilt direction and thus, according to 
the hydrodynamical equations (§ and (§ controls the 
transverse flow of the photon fluid. This can clearly be 


seen in Figs, [^b) and (c) for which the increasing input 
angles lead to increasing flow. The data is consistent 
with Eq. ([^ with An = 0 and flow speeds v^g = 0, 
1.3, 3.0 X 10® m/s [(a), (b) and (c), respectively]. We 
attribute the purely parabolic dispersion to the nonlocal 
nature of the material that, as discussed above reduces 
the effectiveness of the nonlinearity for features smaller 
than the cti,. A nonlocal length for methanol of (Jl ~ 300 
/im has been reported [22] i implying that only modes 
with a much smaller K < 10^ m“^ would appear to have 
a linear dispersion (Fig. [^. 

We underline that in these measurements we are 
imaging the various photon fluid planes inside the actual 
nonlinear medium. This may cause some issues related 
to the fact that we can no longer consider this a linear 
imaging measurement and the nonlinearity may lead 
to deformations of the various object planes [32]. We 
did however, experimentally verify that objects placed 
at the cell input could be imaged through the sample 
and that distortions only appeared to be significant in 
the presence of sharp transitions from very low to very 
high light intensity regions. We therefore expect that 
the images of the weak and smooth oscillations lying 
on top of the more intense laser beam background are 
not significantly distorted and that the overall shape of 
the measured dispersion curve is indeed correct. This 
is corroborated by the fact that the dispersion relation 
Eq. ([^ provides an excellent theoretical fit to the data 
for the zero-velocity flow case and provides a correct 
estimate of the photon fluid flow when the input beam 
is placed at an angle with respect to the imaging axis. 
However, it is clear that a more detailed and precise 
estimate of the low frequency behaviour and the presence 
or absence of superfluid behaviour should be assessed 
by a technique that does not require imaging through 
the nonlinear medium. The results obtained with such a 
technique are discussed in the next section. 

EXPERIMENTAL RESULTS II: PUMP AND 
PROBE TECHNIQUE 

In this experiment, a pump-probe configuration was 
used to imprint a specific wavelength on top of the back¬ 
ground beam as proposed in [1^. We then study the 
wave dynamics for different background intensities. The 
probe beam was extracted from the pump beam using a 
beamsplitter and was attenuated with a half-wave plate 
and a polarising beam splitter. The probe beam is then 
recombined with the pump beam inside the sample: by 
controlling the relative angle of the two beams we create 
an interference pattern of the desired modulation depth 
(« 5%) and relative wavevector K. Both the pump and 
probe beams are loosely focused into the sample using 
a cylindrical lens pair to create an elliptical beam with 
a minor axis diameter of 200 /rm and a Rayleigh range 
longer than the sample (« 20 cm). The interference 
wave pattern (i.e. the photon fluid sound wave) is finally 
imaged at the output facet of the sample and a shift of 
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FIG. 2. a) Experimental layout used for measuring the dispersion relation. A collimated, flat beam is launched through a 
cylindrical sample filled with a methanol/graphene solution and is imaged by a camera/lens system that can be translated along 
the propagation direction of the beam, b) Detail of nonlinear sample showing example beam profiles at different propagation 
distances (i.e. equivalently for different propagation times) as measured in the experiment. An interferometer placed before the 
sample generates a pump and probe beams with a controllable relative angle (i.e. wavelength of the photon fluid excitations). 



K (xlO'’m) 


FIG. 3. Photon fluid dispersion relation with a scanning distance of 1.5 cm, 5 cm, 12 cm (a)-(c). The flow v of the effective 
medium can be controlled by the phase of the background field ( a) Vbg = 0 m/s, b) Vbg = 1.3 x 10® m/s c) Vbg = 3.0 x 10® m/s. 


this pattern along the x direction can be measured as a 
function of the laser power [Fig. |^a)]. The final imag¬ 
ing stage is performed with a 1-to-l imaging f=300 mm 
lens. In the focal plane we place a mask to eliminate 
the four-wave mixed wave, thus suppressing interference 
between counter-propagating Bogoliubov modes and fa¬ 
cilitating the measurement of the phase shift, so as to 
isolate the pump and probe beams. Our pump-probe 
setup is similar to that used in the past to study periodic 
soliton formation [33l [3l] albeit with a lower contrast in 
the pump-probe interference pattern and, as we discuss 
below, our attention is focused on the dynamics of the 
photon-fluid sound waves. 

As the beam propagates through the nonlinear medium, 
the phase velocity of the sound wave is determined by 
vp = n/K with O given by Eq. ([^ and hence is a func¬ 
tion of the nonlinearity An. As the pump beam power 
is increased we therefore expect to observe a spatial shift 
of the sound wave pattern by AS” due to the increase of 
the phase velocity. So by using Eq. ([^, we obtain an 
expression for AS' = AS{An, K): 




|An| 

no 



1 z 


( 8 ) 


By measuring AS for different values of the pump 
power we may therefore use Eq. ^ to estimate both 
An and ap. As a function of the sound wave wavevector 
AT, one can recognize the finite limit AS —>■ An/ng 
determined by the speed of sound in the fluid of light. 
The local vs. non-local nature of the nonlinearity is 
visible for short wavelengths. In the local case, the phase 
shift tends to 0 linearly inl/K with a slope proportional 
to the refractive index change |An|. In the non-local 
case, an extra decay is introduced by the R{K) factor: 
for the case in which which R{K) = 1/(1 -f {(JlK)'^), for 
low K this tends to zero as K~^. As a result, the phase 
shift tends to 0 much more quickly. While a non-zero AS 
is by itself a consequence of interactions, the nonlocal 
effect on the Bogoliubov dispersion is visible in the 
deviation of AS from a linear dependence on 1/K. 


The experimental data of the relative shifts of several 
transverse wavelengths A = 2tt/K between 120 p,m and 
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FIG. 4. (a) Raw data image of the beam at the sample 

output. The sound wave is barely visible due to the low 
(less than 10%) contrast of the amplitude modulation, (b) 
Amplitude profile of the sound wave after substracting the 
pump beam profile for two different input pump powers: low 
power measurement - dotted blue line and high power mea¬ 
surement (showing a shift in the wave, indicated by the ar¬ 
row) - solid red line, (c) Relative shift AS vs. A: the solid 
green line shows the predicted shift for a local nonlinearity of 
An = — 7.6x 10~® - the corresponding healing length is ^ = 85 
fj,m. The solid blue line includes a nonlocal nonlinearity with 
An = —7.6 X 10“® and nonlocal length = 110 fim. The 
shaded gray area (A > 10^) highlights the region in which, 
even in the presence of nonlocality, the medium acts as a su¬ 
perfluid. The black squares indicate the measured shift values 
at various wavelengths, several of which (i.e. for A > 0.7 mm) 
lie in the superfiuid region. 


1.8 mm at a fixed laser ^wer of 28 mW (fluence, 1 
W/cm^) is shown in Fig. 3 (solid black squares). For 
comparison, the predicted shift in a medium with a 


purely local (i? = 1 ) nonlinearity (An = —7.6 x 10“®) 
according to Eq. ([8) is also plotted (solid green line). 
In general, we found that it is not possible to obtain a 
satisfactory agreement with the experimental trend for 
any value of An and a purely local nonlinearity. 

On the other hand, if we take into account the nonlocal 
factor R{K) = 1/(1-|-iT^cr|) in Eq. we obtain a very 
good agreement with the data within = —7.6 x 10“® 
and (Ti = 110 /rm (solid blue line). 

A remarkable feature of these measurements is that for 
the longer wavelengths, i.e. for A > 0.7 mm the observed 
shift appears to saturate and, within an 8% margin, is 
the same for all wavelengths. In other words, the phase 
velocity of the waves, Cg = AS/{Lno/c) = 3.8 x 10® m/s, 
does not depend on wavelength in this region (shaded 
in gray in the figure) and the dispersion relation is 
dominated by the linear behaviour of the Bogoliubov 
term. This is a clear evidence of the linear dispersion 
of collective excitations typical of a superfluid or, in 
the language of general relativity, of Lorentz invariant 
propagation that is crucial for future analogue gravity 
experiments. Naturally, the the observation of superfluid 
features in a system at room-temperature is in itself of 
great interest and future experiments will be directed at 
deepening our understanding and control over the flow 
and wave propagation in this nonlocal photon superfluid. 
A crucial point here is that a standard Z-scan measure¬ 
ment would not be adequate for characterising such a 
nonlocal nonlinearity. Z-scan relies on translating the 
sample through the focus of a lens, i.e. through a region 
in which the transverse length scale of the beam is 
changing. The nonlinearity analysed here shows a clear 
dependence on the transverse beam dimensions through 
the R response function. In a Z-scan measurement the 
material would therefore respond with different effective 
nonlinearities at different measurement positions, result¬ 
ing in a distorted Z-scan trace. We explicitly verified 
that this is indeed the case (data not shown). 

CONCLUSIONS 


We have reported a nonlinear optical experiment 
that aims at characterizing the physics of a photon 
fluid realized in a cavity-less nonlinear medium with a 
nonlocal optical nonlinearity of thermal origin. With 
a technique used in hydrodynamics, we measured 
the parabolic dispersion relation for small amplitude, 
small wavelength noise excitations and highlighted the 
evidence of a hydrodynamic flow in the presence of a 
linear phase term on the transverse beam profile. With 
a pump-probe technique we investigated the larger 
wavelength region by measuring the change in the phase 
velocity of collective excitations as a function of the 
pump intensity. In this way, we estimated the thermal 
nonlinearity to An = —7.6 x 10“® and the nonlocal 
length aL = HO /im. This measurement technique also 
highlighted a superfluid behaviour of the photon fluid, 
which therefore holds promise both for future room- 
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temperature studies of superfluidity and for analogue 
gravity measurements that require a Lorentz invariant 
(i.e. dispersionless) flowing medium. Our measurements 
indicate that there is a readily attainable wavelength re¬ 
gion in which such measurements may be feasible in the 
near future. However, the sound speed in the superfluid 
regime demonstrated is relatively low Cg = 3.8 x 10® m/s 
and will require very long samples in order to observe 
significant wave propagation. Future work will therefore 
be directed at increasing the sound speed, e.g. by use of 
materials with higher nonlinearities. On the other hand, 
the degree of nonlocality can also be tuned by modifying 
the sample geometry [23] or absorption [35] and this 
opens opportunities for performing related experiments 
over a broader wavelength range. A crucial and natural 
demonstration of superfluidity will involve the study of 
the transition from laminal flow arond an obstacle at 
subsonic flow speeds to turbulent dissipative flow with 
the formation of quantised vortices at supersonic flow 
speeds [iniiMi- Furthermore, such media show promise 
for studies of hydrodynamical turbulence in the presence 
of nonlocal interactions and the study of collective 


dynamics of incoherent waves. Finally, extension of our 
experimental apparatus beyond the CW measurements 
shown here so as to measure also time-dependent 
fluctuations could provide access to quantum fluctua¬ 
tion features [H] such as analogue Hawking radiation |3]. 
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